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The staggered fermion determinant is complex when the quark chemical potential fi is nonzero. 
Its fourth root, used in simulations with dynamical fermions, will have phase ambiguities that 
become acute when Re fi is sufficiently large. We show how to resolve these ambiguities, but our 
prescription only works very close to the continuum limit. We argue that this regime is far from 
current capabilities. Other procedures require being even closer to the continuum limit, or fail 
altogether, because of unphysical discontinuities in the measure. At zero temperature the breakdown 
is expected when Re fi > where m, is the pion mass. Estimates of the location of the 

breakdown at nonzero temperature are less certain. 



The direct study of quantum chromodynamics (QCD) 
at nonzero chemical potential \i faces the considerable 
hurdle presented by the complex measure in the euclidean 
path integral (for a review, see Ref. lj). This requires the 
adoption of various stratagems in order to make Monte 
Carlo analysis possible. These stratagems include Tay- 
lor expansions around \i = 0]; analytic continuation 
from the imaginary fi axis (where the measure is real 
and positive) |jfl; and multiparameter reweighting using 
real measures that push the complex phases into the ob- 
servables 

II H IS IE 1 All these methods envisage a 
well-defined path integral with complex measure at real, 
nonzero values of fj,. While the first two methods never 
actually calculate fermion determinants for Re fi ^ 0, the 
third method does calculate ratios of these determinants 
as configurations are analyzed. For a different reweight- 
ing method based on the canonical ensemble, see Ref. 0. 

Calculations at fi ^ are often performed with 
staggered fermions. These are contained in a single- 
component fermion field that makes for efficient calcu- 
lation of the fermion determinant and Green functions; a 
further advantage is a chiral symmetry that is only bro- 
ken softly by the mass parameter m. A disadvantage of 
staggered fermions is residual doubling, wherein a sin- 
gle lattice field actually represents four physical flavors 
(nowadays called tastes) that become degenerate in the 
continuum limit. All these tastes affect the dynamics of 
the theory, as is easily seen by considering fermion loops 
in a diagrammatic expansion of the determinant. In or- 
der to represent a single flavor of fermion with mass m, 
then, one takes the fourth root of the determinant. Thus 
one hopes to arrive at two light flavors and one strange 
quark in the continuum limit. 
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The QCD partition function on a lattice, with two light 
quarks plus the strange quark, is given by 

Z = J *VUe~ s ^ u) det[C(C/)+m i ] 1 / 2 

xdet[D(U) + m s ] 1/4 , (1) 

where D(U) is the staggered- fermion hopping matrix, in- 
cluding factors of exp(±/i) on the timelikc links; mi is 
the common mass of the two light quarks and m s is the 
strange quark mass. When Re/i = 0, the matrix D(U) 
is anti-hermitian. The eigenvalues of D + m then come 
in conjugate pairs, rj — m ± iE, so the determinant is 
real and strictly positive (for m ^ 0). If we take the 
fourth root to be real and positive as well, the measure 
in Z is a continuous functional of the gauge field U. We 
will assume here that taki ng t he fourth root is a valid 
procedure for \i = [13, IH 111 • 

When Re [j, =^ 0, on the other hand, the spectrum of 
D + m spreads out into the complex plane and the de- 
terminant becomes complex. Taking the fourth root of 
the complex determinant introduces phase ambiguities. 
A prescription must be found for resolving these ambi- 
guities. The problem first appears when fi > m, when 
rare configurations admit eigenvalues near zero. At some 
larger value of //, near-zero eigenvalues become common 
in typical configurations. Analysis of the e-regime sug- 
gests that in the zero-temperature theory this happens 
rather abruptly when fi w |13| : at T > this point 

is believed to move towards higher \i Q . 

We present below a sensible prescription for taking the 
fourth root when Re /i ^ 0. This prescription requires 
that we work near the continuum limit. The breakdown 
of the method on coarser lattices, in fact for any realistic 
value of the lattice spacing, is quite general and indicative 
as well of the failure of alternative rules for dealing with 
the fourth root. The discussion carries over to the square 
root used for the two light quarks. 
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An unambiguous definition of the phase of the fourth 
root requires understanding taste symmetry and its 
breaking at finite lattice spacing. To this end, let us de- 
scribe the IR and UV regimes of the complex spectrum 
of D + m near the continuum limit (see Fig. QJ. The 

, 




FIG. 1: Sketch of complex eigenvalues rji of D + m when Re/i 
is somewhat larger than The dashed line is Re 77 = m. 

Eigenvalues closest to the real axis are clustered into quartets 
with splittings of 0(aA 2 ); typical spacing between quartets is 
oc V -1 / 2 (see text). Farther from the real axis the quartet 
structure disappears. The eigenvalues come in pairs (77, 2m — 

eigenvalues rji that fall in the IR regime will form ap- 
proximately degenerate, well distinguished quartets, as 
dictated by taste symmetry. These will include eigen- 
values near the QCD scale, \r)i\ « A. For UV eigen- 
values, with |7/j| — > o , large taste violations are un- 
avoidable. Since, however, both m and /_t are now small 
compared to these eigenvalues are dominated by the 
anti-hermitian D(/i — 0), and so their phases should be 
near ±7r/2, with deviations that are 0(am,a[x). We ex- 
pect that there will be a large overlap region where the 
eigenvalues rji satisfy both conditions: they should fall 
into approximately degenerate quartets and their phases 
will be near ±7r/2. 

Now the prescription for the fourth root of the determi- 
nant is straightforward. For the IR quartets, we proceed 
one quartet at a time, and we take the fourth root of their 
product to have a phase angle equal to the phase of the 



center of the quartet. Thus we construct a single- flavor 
theory, with one eigenvalue in place of each quartet. For 
the UV eigenvalues, 1 there is no need for subtlety; since 
the eigenvalues lie near the imaginary axis, we can take 
the arguments of their fourth roots to lie near ±7r/8. 

The fourth root so constructed will usually be a well- 
defined, continuous functional of the gauge field, with a 
slowly varying phase angle. An exception will arise if 
an eigenvalue quartet comes too close to the origin. If, 
in fact, the four eigenvalues frame the origin, the aver- 
age phase of the quartet becomes meaningless and any 
definition will fluctuate rapidly as the eigenvalues move. 
Such an occurrence will, however, be rare. The typical 
splitting 8 within a quartet will be 2 0(aA 2 ), and so the 
probability of finding an eigenvalue this close to zero is 
proportional to a 2 . A more detailed estimate follows. 

First, let us discuss orders of limits. At fJ, = 0, the 
relative taste splittings inside a quartet, \rji — rij\/\rj\, are 
bounded from above by aA 2 /m. If we want the taste 
splittings of the pseudo-Goldstone bosons to be small, we 
must take the continuum limit before the chiral (m — > 0) 
limit 0. Let us now assume Re/j ~ A. The eigenvalue 
density near the origin in the complex r\ plane scales like 
VA 2 with the volume V of the lattice. Thus m gets re- 
placed by the 0[(VV A) -1 ] mean distance between quar- 
tets (the same estimate applies to the typical magnitude 
of the eigenvalue nearest to the origin), and \rji — T)j\/\fj\ 
comes out to be 0(a\/V A 3 ) for the near-zero quartets. 
Now it is the product aVV A 3 that must remain small 
for the sake of good taste. If not, the splittings within a 
quartet will be as large as the splittings between different 
quartets, and our prescription will break down. When 
approaching the continuum limit, the volume must be 
limited such that ay/V A 3 — > as well. 

With the eigenvalue density estimated as above, the 
fraction of quenched configurations containing an eigen- 
value quartet that frames the origin is 0(a 2 VA 6 ). This, 
then, is the order of the error introduced by either tak- 
ing or dropping these configurations. Reweighting by the 
determinant provides further suppression by 0(ay/V A 3 ); 
this is the ratio of fourth roots of the said eigenvalue 
quartet and of a typical near-zero quartet. The system- 
atic error thus becomes 0[(aVV A 3 ) 3 ]. 

Is it really necessary to extract individual eigenvalues 
and to group them into quartets? The procedure appears 
cumbersome, but other prescriptions that come to mind 



1 Renormalization-group analysis llll suggests that the log of the 
UV part of the determinant contributes only a local correction to 
the effective action of the gauge field; it can even be discarded, 
if one compensates by renormalizing the gauge coupling. 

2 This can be seen from the representation of the staggered op- 
erator in taste basis, where the taste breaking comes from a 
dimension-5 operator that shifts eigenvalues within a quartet as 
r? — » 77±<5, where S is fixed and 0(a). This is consistent with the 
usual picture of 0(a 2 ) discretization errors because the product 
of eigenvalues will shift by 0(5 2 ). 
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produce larger errors close to the continuum limit, or 
even fail altogether. If one were to calculate the determi- 
nant T> in full and only then to take its fourth root A in 
the cut plane |Arg2?| < ir, the result would be restricted 
to the wedge |Arg A| < ir/4. This is inconsistent with an 
ordinary single-flavor theory, whose determinant is free 
to wander the entire complex plane. 

If one were to take the fourth root of each individual r\i 
in the cut plane before multiplying them together, there 
would be frequent jumps of ir/2 in the phase of A as single 
eigenvalues cross the cut. Here the systematic error can 
be estimated: It would be 0(aVA 5 ), which is the fraction 
of configurations in which a quartet is bisected by the 
negative real axis, close to the continuum limit. [The 
probability that none of the 0(VA 4 ) IR eigenvalues fall 
in a strip of width aA 2 along the negative real axis scales 
like (1 - aA) VA ps 1 - aVA 5 .] This estimate is only valid 
if one has verified the existence of well-defined quartets in 
all configurations. Moreover, this error is can be fatally 
large even if the figure of merit ay/V A 3 is small. 

Are current calculations anywhere near to realizing the 
quartet structure of IR eigenvalues that will make calcu- 
lation of the fourth root consistent? We suspect that the 
answer is negative. Simulations at fi = show nice taste 
degeneracy in the lowest eigenvalues only when the lat- 
tice cutoff is a -1 ~ 2 GeV or larger, and only if highly 
improved gauge and staggered actions are used 0. The 
state of affairs is much less satisfactory in existing sim- 
ulations with Re/i > 0, where the cutoff is at 1 GeV 
or less. Even in the fj, = simulations, significant taste 
violations will reappear if one looks beyond the lowest 
few quartets. While this does not hamper simulations at 
/i = 0, for the /x > case we do require a distinct quartet 
structure all the way up to \rj\ ~ A and beyond. 

In general terms, what would it take to define a fourth 
root that is single-valued and continuous? Consider a 
closed trajectory in gauge-field space that avoids the ze- 
ros of the fermion determinant D. When this trajectory 
is mapped into the complex T> plane, the winding number 
around the origin must be a multiple of 4. In practice, 
this will never be guaranteed, precisely because of taste- 
symmetry breaking. We believe that current simulations 
are plagued by eigenvalues with no particular multiplet 
structure. Problems will arise when any single eigenvalue 
circles the origin. Depending on the treatment given the 
fourth root, either the eigenvalue will cross a cut, caus- 
ing a jump in the fourth root; or the eigenvalue will move 
onto another Riemann sheet, giving a multi- valued root. 

[We note that, for m = 0, two eigenvalues will always 
circle the origin simultaneously, because they appear in 



(17, —rj) pairs in that case (see Fig. QJ. It would then be 
possible to define a square root by requiring continuity, 
because the winding number of the determinant around 
the origin would always be even. This only works, how- 
ever, for m = 0, and it will not help for the fourth root.] 
We illustrate this problem in two dimensions. We shall 
exhibit eigenvalues near the origin for selected gauge 
fields and show how small changes in the gauge field 
can take a single eigenvalue around the origin and thus 
through any cut that can be defined. The staggered 
fermion action on a lattice of L t x L x sites is 



S 



4>{D + m)ip 



t=l x=l 



+2rmp tx ij) tx 



(2) 



A change of variables [Tg has removed the chemical po- 
tential from all links of the lattice except for those that 
connect t — L t to t = 1. First we choose a gauge field 
that is zero except for a constant Polyakov loop U on the 
last time-link. This is equivalent to writing the free oper- 
ator with the modified antiperiodic boundary conditions 



A Fourier transformation from x to p gives 
Lt { 1 

s = } 2 i^pipt+hp ~ ipt+i,pipt, P ) 

t=i v 



+ [m+ (-l)*ijt>] ipt,pipt,p 

t,v p 



(4) 



where p — sin p. The momentum takes the values 
p n = (2n+ \ )tt/L x . We must keep careful track of multi- 
plicities. If we take L x to be a multiple of four, then none 
of the momenta p = 0, 7r/2, tt, —tt/2 are in the spectrum. 
This assures that each possible value of sin 2 p has a mul- 
tiplicity of four. The fermion matrix in the momentum 
basis is given explicitly by 
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M(p) = 



m — ip 
_ 1 
2 
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We have chosen £7 diagonal, whereupon the determinant 
factors into N c color factors, and selected one of the col- 
ors. 

M is not anti-hermitian; it satisfies instead 
M'(m,n) — —M(—m, —/!*). The right-eigenvectors of 
M may be found via an alternating ansatz, 



in 



ae 



ial 



f3e i0t . 

We find the eigenvalues 



t odd, 
t even. 



rj± = m ± i£, 



(5) 



(6) 



where £ = ^Jp 2 + lj 2 , with u> = sin a;. The boundary 
conditions lead to the quantization of Re0, viz., 



Re 0„ = Im \x 
and further require 



(2n + 1)tt- (f> 



(7) 



Im0„ = —Re fi. 



The eigenvalues depend on p only via p 2 = sin 2 p, and 
thus the complete set of eigenvalues shows a four-fold 
degeneracy. 



When Re^i 



> 



[more precisely, cosh(2Re^) > 1 



2m ] , the eigenvalue ry + can be zero if £ is pure imaginary, 
i.e., £ = im. This can be accomplished by giving the 
Polyakov loop U a value corresponding to 



L t Im /i + 7r (mod 2n) ; 



(9) 



the zero eigenvalue is then found by choosing p appro- 
priately and taking n — so that Rew = 0. A small 
change in U from this value gives us a four-fold degener- 
ate eigenvalue near but not at zero; in fact, p cannot be 
chosen to give an exact zero in the first place unless L x 
is infinite. 

We can add a gauge field to break the degeneracy by 
pushing the four eigenvalues in different directions. Con- 
sider a spatial gauge field that, for all x, takes values 



JA 



and Ut 



for odd and even t respec- 



tively. (We assume that these spatial links commute with 
the constant Polyakov loop.) In Eq. (0J, this replaces 
p — smp by p = sin(p + A a ) and by p e = sin(p + A e ) for 



lmT| / 






ReT| 



FIG. 2: Breaking of the fourfold degeneracy and transport 
of one eigenvalue around the origin. The small filled circle 
marks the original (fourfold) eigenvalue r\. The two upper 
trajectories have been slightly displaced for clarity. 



odd and even t respectively. The spectrum can again be 
calculated exactly. The new eigenvalues are 



r)'±=m+ i(p e 



Po )/2±i£\ (10) 

j 2 . The conditions (JJJ) and 
on lj are unchanged. If we expand to first order in 
A a and A e , we obtain rj'^ = rj± + 5i}±, where 



(g) where £' = sj \p Q + p e ) 2 / '4 + 



6r]± 



i cosp 



(A e -A )± S -^(A 



A e 



(11) 



Equation contains both cosp and sin p. This means 
that the four-fold degeneracy among p, — p, n + p and 
7T — p is completely lifted by the gauge field (see Fig.GJ. 3 

To see how we can make a single eigenvalue circle the 
origin as m Fig. let us be gin again from the degenerate 
case. The eigenvalue quartet closest to zero will be of 
order l/L x -c 1. We turn on A Q , then A e , and then 
turn them off in the same order. The closed trajectory 
of only one eigenvalue circles the origin once; hence also 
the determinant circles the origin once. 



3 When the gauge field is zero the fourfold degeneracy is a result 
of taste symmetry plus charge conjugation. The constant gauge 
field A e = A breaks the latter, leaving the degeneracy between 
p and p + 7r. Only when A e ^ A is the taste symmetry broken, 
as can be shown with the method of Ref. llSl . 
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A local fluctuation in the gauge field can accomplish 
the same lifting of the degeneracy and transport of eigen- 
values. For instance, beginning with the degenerate 
case, choose a single spacelike link at (t, x) and set 
U} x = 1 + iA, with A chosen to be diagonal in color. 
The first-order change in a zero eigenvalue is 

tAcosp f-l-^sinp, todd, 
= j j x \ ■ ( 12 ) 

[ 1 ~Wi sin P> * even- 
Again, the result depends on both cosp and sinp, so that 
the four-fold degeneracy is lifted. 

Our study of two-dimensional QCD exhibits nondegen- 
erate eigenvalues and shows how they wind around the 
origin. The fourth root of the determinant will be prob- 
lematic here; the same will apply to any representation 
of the determinant as a product, 

det(U + m) = (13) 

i 

where the factors are complex functions of the gauge field, 
if the fourth root of the determinant is taken by multi- 
plying the fourth roots of the individual £j's. 

An example is the formula used by Fodor and Katz 
jHQ) where £j = e ifAI — \f, the \ are eigenvalues of a 
non-hermitian matrix on a three-dimensional sublattice. 
In our two-dimensional example, again one finds that by 
turning on the gauge field A ey0 the degeneracy of the A; 
is completely lifted. 

The method of Ref. [f| guarantees continuity of the 
fourth root of each as a function of fi, starting from a 
point where Re \i = 0. This is not the same as continuity 
as a function of the gauge field, which, as we have seen, 
is the central issue. Each £j carries within it a tangle of 
information coming from all length scales; we do not even 
know if the £j's will fall into well-distinguished quartets 
close to the continuum limit. If they do, it would follow as 
above that the systematic error of this prescription is due 
to £i's located a distance of order a from the negative real 
axis. As noted, this gives a systematic error of 0(aVA 5 ). 
If there are no quartets then the error will be larger. 

We can make our estimates semi-quantitative for the 
simulation parameters of Ref. |3| . It is convenient to iden- 
tify the QCD scale A with the critical temperature T c at 
/i = 0. Since the temporal size of the lattices used was 
in the neighborhood of 1/T C , in these units the lattice 
spacing was a ~ 1/L t = 1/4; the largest volume used 
was (12/4) 3 = 27. As we have seen, if the eigenvalue 
spectrum is to exhibit a taste-quartet structure then the 
figure of merit ay/V A 3 must be small. For these param- 
eters, however, we arrive at a value of 1.3. This returns 
us to a picture of no quartet structure and hence an in- 
tractable error due to phase discontinuities. 

The remaining question regarding Ref. 8] is whether 
the values of \i studied put the theory into the danger- 
ous regime where the distribution of eigenvalues rji or £j 



reaches the origin. In fact Re ll was varied up to and be- 
yond the putative critical end point at \i c ~ 120 MeV. 4 
This is well above the value of 771^/2 for the realistic 
light-quark masses used, but, as we have noted, the crit- 
ical value of n is probably larger when T > 0. The fact 
that the pseudo-critical temperature T c (ji) determined in 
Ref. agrees over a substantial interval with that cal- 
culated via analytic continuation |3J is indicative of such 
a shift, perhaps as far as [i ~ 100 MeV. Direct study of 
the spectral density would settle this question. 

Let us summarize our findings. While there are new 
complications, taking the fourth root when Re \i =/= 
could in principle be on a footing similar to that of the 
fj, = case. A requirement is that one work close enough 
to the continuum limit and set the volume such that taste 
multiplets have reappeared in the eigenvalue spectrum 
all the way up to \t)\ ~ A and beyond. If one tracks the 
quartets and assigns their phases carefully, then config- 
urations with quartets that frame the origin will bring 
a systematic error of 0[(aWA 3 ) 3 }. If one calculates 
fourth roots of individual eigenvalues and ignores their 
quartet structure, then quartets that are bisected by the 
cut cause an error of 0(aVA 5 ). In practice, however, 
the demonstration of the necessary quartet structure ap- 
pears to be far beyond any realistic simulation even at 
Re ll = 0, let alone when Re fi ^ 0. 

In simulations done at purely imaginary \i |3j the valid- 
ity of the fourth root is on the same footing as at ll = 0, 
provided the continuum extrapolation is done before at- 
tempting the analytic continuation in /i. Otherwise, on 
a finite lattice, the integration over the compact gauge 
field commutes with the analytic continuation in ll. The 
method would then be tantamount to some prescription 
for taking the fourth root of the determinant on the in- 
dividual configurations. The nature of this prescription, 
and its relation to the ones compared here, is beyond the 
scope of this paper. 
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4 This is the quark chemical potential; Ref. 0] gives the baryon 
chemical potential, fig = 3^i. 
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